Weak non-linear surface charging effects in electrolytic films 



D.S. Dean(i) and R.R. Horgan^^) 
(1) IRSAMC, Laboratoire de Physique Quantique, Universite Paul Sabatier, 

118 route de Narbonne, 31062 Toulouse Cedex 04, France 
(2) DAMTP, CMS, University of Cambridge, Cambridge, CB3 OWA, UK 
E-Mail:dean@irsamc.ups-tlse.fr, rrh@damtp.cam.ac.uk 
(Dated: 7 October 2002) 

A simple model of soap films with nonionic surfactants stabilized by added electrolyte is studied. 
The model exhibits charge regularization due to the incorporation of a physical mechanism respon- 
sible for the formation of a surface charge. We use a Gaussian field theory in the film but the full 
non-linear surface terms which are then treated at a one-loop level by calculating the mean-field 
Poisson-Boltzmann solution and then the fluctuations about this solution. We carefully analyze the 
renormalization of the theory and apply it to a triple layer model for a thin film with Stern layer 
of thickness h. For this model we give expressions for the surface charge (j{L) and the disjoining 
pressure Pd{L) and show their dependence on the parameters. The influence of image charges nat- 
urally arise in the formalism and we show that predictions depend strongly on h because of their 
effects. In particular, we show that the surface charge vanishes as the film thickness L — > 0. The 
fluctuation terms about this class of theories exhibit a Casimir-like attraction across the film and 
although this attraction is well known to be negligible compared with the mean-field component for 
thick films in the presence of electrolyte, in the model studied here these fluctuations also affect the 
surface charge regulation leading to a fluctuation component in the disjoining pressure which has 
the same behavior as the mean-field component even for large film thickness. 



I. INTRODUCTION 

Many situations in colloid, polymer and inter-facial science involve charged objects interacting in electrolyte solu- 
tions. In the case of two interacting membranes or soap films one encounters the electric double-layer. The mean-field 
theory of such experimental configurations is the Poisson-Boltzmann theory |^ for which a surface boundary con- 
dition must be determined. For example, the surface charge or surface potential can be given, or a relation between 
the surface charge and surface potential may be specified. Such a relationship arises in a charge regularized model 
where the surface charging mechanism is derived from a microscopic description of the chemistry and geometry at 
the interface. Charge regularized models have the property that, owing to the thermodynamic nature of the charging 
mechanism, the surface charge changes as the distance between the two surfaces varies. The variation of the surface 
charge with the inter-surface distance will also change the effective interaction between the surfaces and consequently 
also the disjoining pressure. Within the film, the mean-field Poisson-Boltzmann theory is only sensitive to the elec- 
trostatic properties of the electrolyte; the chemistry and effective sizes of the ions in the system only enter into the 
description of the surface charging process. For example experiments ||^ shows that the disjoining pressure increases 
with increasing hydration radius of the counter-ions in ionic soap films. This effect can be explained by the fact that 
as the counter ion radius increases the capacity for it to approach the surface and screen the surface charge is reduced. 

To take into account bulk surface tensions of electrolytes one must resort to a microscopic description of the physics 
at the interfaces. Depending on the nature of interface one may have specific adsorption due to chemical effects 
which can be taken into account via the law of mass action |Q, ||, ^, ^, ^ or at a more statistical mechanical level by 
introducing an external potential at the surface which models the chemical liaison involved 0. In addition, 

there are other forces which come into play which are not present in the standard Poisson-Boltzmann approach. 
These are effects due to fluctuations in the electromagnetic field which can be identified with van der Waals forces 
0, [13[ g. In bulk, the non-zero frequency van der Waals forces have little effect on the static ionic distributions 
since the relevant frequencies are too large and their contributions can be decoupled from those of zero frequency. 
However, it was pointed out in that near interfaces the non-zero frequency van der Waals or dispersion forces 
can be important and depend strongly on the polarizability of the ions involved and hence are ion specific. The 
zero frequency van der Waals forces, which correspond to thermal fluctuations in the electrostatic field, do strongly 
influence ionic distributions and so do modify the surface charge. 

When there are spatial variations in the dielectric constant image charges arise fl^ . In the field theoretic approach 
adopted in this paper image charges and their effects are naturally and systematically included by taking into account 
the fluctuations of the electrostatic field. 

In section || we describe the field theory model we use to study the effects of non-linear terms on the surface charging 



mechanism whilst retaining the free field theory description for the bulk electrostatic fields. This corresponds to using 
linear Debye theory in the bulk with fugacity /i but with the fully interacting description of the sources for the charging 
mechanism of the surfaces. The model is applied to a triple layer thin film as an idealized model of the surfaces where 
there is a Stern layer of thickness h from which all ions are excluded. In section III we discuss the mean-field solution 
to the field theory for the thin film for which the non-linear interactions at the surface determine the source terms 
in the mean-field equation. In section |^ we give a detailed description of the effect of field fluctuations about the 
mean- field solution using the Schrodingcr kernel approach developed in an earlier paper psf . We give predictions for 
the surface charge cr(L) and the disjoining pressure Pd{L) as a function of film thickness L and show under reasonable 
assumptions that cr{L) ^ as L ^ and that the mean-field prediction for the large L tail of Pd{L) is modified by 
one-loop corrections, both decaying like exp(— mL) for large L, where m is the Debye mass . In section ^ we present 
a number of example graphs and discuss their salient features. In particular, it is clear that the effects of image 
charges, which arise naturally and systematically in our formalism, are very strong and that qualitative predictions 
depend sensitively on the value of /i, the thickness of the Stern layer. This demonstrates that a realistic model for the 
structure of the surface and the Stern layer in particular is necessary for a quantitative study. We also present some 
conclusions in this section. 



II. MODEL 



We consider an idealized model of a thin film made with nonionic soap adapted to the experimental set up used to 
measure the disjoining pressure as a function of the film thickness |lo[ The film shown Fig. (|l|) consists 

of two parallel surfaces of area A the interior of which is filled with a monovalent electrolyte solution such as NaCl 
in water with bulk dielectric constant e. The exterior of the film is a dielectric medium of dielectric constant eg, for 
example, air. The perpendicular distance between the two surfaces is denoted by L -I- 2h. The region of thickness h 
is the Stern layer from which the largest ions in the system are excluded and h can be taken to be the radius of the 
largest ion type in the system, which is here chosen to be the anion (in most physical systems it is the cation with the 
largest radius due to hydration). If the radius of the cation in solution is h' and h' < h then the cation is excluded 
from a region of width h' from the surface but can be present in the region [—{h — h'),0] where there are no anions, 
thus leading to an effective surface charge in that region. Strictly speaking, just outside the film is the surfactant 
layer which in general will have a different dielectric constant to that of the exterior and the aqueous interior. Here 
for simplicity the presence of the nonionic surfactant is neglected. A version of this model with a single surface plus 
bulk is commonly used to model the surface properties of electrolyte solutions 22 . In experiments the 



thickness L of the film may be varied by applying an external pressure in the cell containing the film and its bulk. 
The grand partition function for this model system may be expressed as a functional integral 



s= / d[<p]cxp{sm (1) 

where 



S[(j)] = / /3e(x)(V0)2dx + 2n [ cos {(3e(j)) dx 



+ /i+ / {S{z) + 6{L ~ z))exp{i(3e(j)) dx , (2) 

JlxA 

where e is the electron charge, A is the area of the film and f3 — l/kgT. The fugacities of the anions and cations 
are taken to be equal and denoted by fi as are their densities denoted by p. The total length in the z direction 
perpendicular to the film surfaces is T + L + 2h where T denotes the total length external to the film. In the region 
z € [— /i, L-\-h] the dielectric constant, which is a function only of z, is given by e{z) = e and outside the film e(z) = eq. 

The above field theoretical formulation can be obtained directly from Quantum Electrodynamics |2^ by retaining 
just the electrostatic potential field in the QED Lagrangian and coupling it to the distribution of ion charges. Alterna- 
tively, it can be obtained by standard field theory techniques based on the Hubbard-Stratonovich transformation for 
a monovalent Coulomb gas. The formulation takes into account both the Coulomb interactions between the ions and 
the zero frequency van der Waals forces due to fluctuations in the electrostatic potential Q. These zero frequency 
van der Waals forces are particularly relevant when there are variations in the dielectric constants of the system since 
they naturally and systematically include the effects of image charges. This is especially the case for aqueous soap 
films in air where e/eo ~ 80. The density operators for the cations and anions can be shown to be 

P± (x) = exp (±ie/30(x)) (3) 
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The last term in Eq. (g) represents a highly localized affinity for the cations to be at the surfaces z = and z = L oi 
the film, and is responsible for the generation of a surface charge. There are various mechanisms leading to affinities 
for ionic species at interfaces, ranging from chemical affinity to steric and entropic effects. Here this term arises in 
the following approximation. Since there can be cations in the region [—{h — h'),0] and the corresponding region 
[L, L + {h — h')] there is, in addition to the first two terms of the action Eq. (0), a surface term 



Yi = fi exp (iPecj)) dx + fi exp {iPe(f>) dx . (4) 

J[-{h-h'),0]xA J[L,L+{h-h')]xA 

This is the integral of the density operator given by Eq. (^) for the cations over the regions [—{h — h'),0] and 
[L, L + {h — h')]. These are the regions in the Stern layers which may be occupied by the cations but from which the 
anions are excluded, li h — h' <^ lo, where lo is the Debye length which characterizes the scale of the electrostatic 
interactions, then we can approximate S by 

S = jjiih — h') I 5{z) exp {iPecf)) dx + — h') I 5{z — L) exp (iPecf)) dx (5) 
JlxA JlxA 

Comparing with the formula Eq. (||) we see that — ^i[h — h'). In general actions of the type Eq. ^ can be used 
to describe various surface charging mechanisms with the proviso that the region where the surface charge is localized 
has a width much smaller than the Debye length, Id- Whilst this description acts as a motivation for the model it is 
clear that in general there will in general be specific adsorption for the different species near the air water interface. 
To model a specific adsorption for the cations in the surface region [—{h — h'), 0] we take a surface cation fugacity /i^ 
in this region which is greater than the bulk cation fugacity /i. On shrinking the surface region to zero one would then 
have = fj,s{h — h'). There is some evidence for specific absorption of certain ionic species at seemingly chemically 
neutral interfaces; the famous and controversial Jones-Ray dip |24|| in the surface tension of weak electrolyte solutions 
(with interfaces with air) can be explained by invoking a specific adsorption of anions at the interface |25[ , although the 
basic surface exclusion model introduced in |po| and used here cannot explain negative excess surface tensions. Recent 
experimental evidence points towards a specific adsorption of hydroxide ions at air/water and oil/ water interfaces 
p2| . In what follows we treat ^*^_ as an independent variable to account for the more complex charging mechanisms 
which can occur at the interface. 

It is important to note that the fugacities /i and /x* are determined by the values of the physical bulk density p and 
surface density for a single planar surface p^. We then have 

p = /i(cos/3e0) , p*^ = p*^{ex\i{i(3e(j))) , 

where the brackets stand for averaging over the bulk partition function. In mean-field theory this implies p = p, P+ = 
p*^ but this is not true in general and the relationship must be calculated taking field interactions into account. 

In the field theory with the action of Eq. (H) was analyzed in the weak-coupling or Debye-Hiickel limit which is a 
Gaussian approximation where the action is expanded to second order in the field (f). This amounts to the assumption 
that the mean-field densities of cations and anions throughout the film are small enough so that 8ttp{z)1^ < 1, 
where Ib is the Bjerrum length Is — e^/3/47re. Another approach is to solve the full non-linear mean-field equations 
and then calculate the one loop correction which gives the effect of field fluctuations about this mean-field solution 



|g7, 28 . In |g7| the resulting mean-field solution in the case of fixed surface charge was ingeniously expressed in terms 
of special functions which allowed an analytic calculation of the one loop correction. In our case, however, the fact 
that the surface charge must be computed self consistently leads to additional complications and we take a different 
approach. We assume that the Gaussian approximation is valid inside the film but not at the surface where, because 
of the increased charge density, the electrostatic interactions will be stronger and we retain the full non-linear surface 
operators. The resulting theory can be analyzed as before but with an effective L— dependent surface source for 
the field. The theory is accurate to the same order in perturbation theory as before but, in addition, now includes 
non-perturbative surface effects. The corresponding Debye-Hiickel action with non-linear surface terms is 

S'*[0] ^ - \ I I3e{z){\7(t>fdx~-]- [ (3em^(j? dx 

2 J {T+L+2h)xA 2 Ji^yA 



P*+ / iS{z)+6{L- z))exp{i(3e(t))dx + 2pAL , (6) 



LxA 



where m = y/SnplB = is the bulk Debye mass, is the Debye length and Ib is the Bjerrum length defined 
above. The weak coupling limit corresponds to ml b < 1, and the Gaussian approximation will be valid throughout 
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the film so long as the local or effective mean-field Debye mass m(z) = ^y8^Tp{z) Is does not become so large that this 
weak coupling condition is violated. This condition is 8TTp{z)lg < 1 as stated earlier. We note that since p(z) > pbuik 
this necessarily requires that Snpbuik < 1 ■ 

Proceeding with the approximation scheme described above we have the expression 

d[0]exp(^*[0]) (7) 

for the grand partition function. The mean-field Poisson-Boltzmann equation is obtained from the saddle point of the 
action S*: 

. = (8) 



<5<^(x) 



with p*^_ replaced by p*^. This is the correct procedure since perturbative corrections due to the interactions of field 
fluctuations about this mean-field, or tree level, solution will then relate p^ and p"^ self-consistently order by order. 
In principle, the same approach applies to the bulk quantities p and p but since we are limiting the current analysis 
to a Gaussian theory in the bulk, the distinction need not be made at this stage. 

Taking into account the Gaussian fluctuations about the mean-field solution gives the perturbative correction to 
one-loop for which we have 

S « exp {S*[^e]) J exp Q J dxdy^^^^^U>'(x)0'(y)^ , (9) 

with, again, p'^ replaced by p^. All other terms are treated as interactions to be analyzed by perturbation theory. 

The grand-potential per unit area of film J can be separated into a mean-field contribution plus the zero frequency 
van der Waals contribution coming from the fluctuations. We write 

J = j*^^ + r""^ (10) 

where 

(11) 

and 



III. MEAN-FIELD THEORY 

The mean-field equation is obtained as usual by looking for an imaginary solution to Eq. (^) <j)c = iip where if) is 
real and corresponds to the mean- field electrostatic potential The resulting equation for is 

/3V-eVi/;-m2/3e^ + p;/3e((5(z)+(5(L-z))exp(-/3ei/;) =0, (13) 

within the film and outside the film one has 

/3V • eoV?/' = . (14) 

This mean-field equation has the form of a standard linearized Poisson-Boltzmann equation but with non-linear 
boundary terms with p*^ is replaced by its mean-field approximation p*^_ as explained in the previous section. 

The solution for ■;/; is by symmetry only dependent on z and symmetric about the mid-plane of the film aX z — L /2 
and so we choose the solution i^{z) — C(L)cosh (m(z — -j)) inside the film. Outside the film Eq. ( p^ gives that 
dip/dz = 0, which is the condition of electro-neutrality of the mean- field solution within the film. Integrating the 
mean- field equation between z = 0~ and 2 = 0"*", and using the condition of electro-neutrality, one finds 

e2u=o+ = -ep*+ exp (-/3e^(0)) (15) 
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Defining D{L) = C{L) /3e cosh(m_L/2) gives the non-linear self consistent equation determining D to be 



D{L) = a cotli(^) eM~D{L)) (16) 

where a = mp\l2p is dimensionless. In the Stern layer model a = m{h—h')/2 = {h — h')/2lD- As stated in Section || 
it is when Id 3> {h — h') that the formulation in terms of a surface charge is valid. In the study of this particular model 
one has that a <C 1, which also implies small surface charges compatible with the use of the quadratic approximation 
within the film. Then Eq. (|l^) may be formally solved as a power series in a coth(^^) using standard series inversion 
techniques from complex analysis. We find that 

DiD^a coth f ^) V (-1)" fa coth (^4-\X (!i±i^ (17) 



V 2 / 'V V 2 

We note that in general surface charge regulated models | ^ , even if there is an exact solution to the Poisson-Boltzmann 
equation (or its linearized form) in the bulk, one must determine the surface potential via a transcendental equation 
relating the surface charge to the surface potential which in general may be solved numerically by iteration ^ or by 
linearizing the boundary equation |^ . Fortunately in the case studied here we have an explicit series solution to the 
boundary equation. The disjoining pressure of a film Pd{L) is the difference between the film and bulk pressures. In 
the grand canonical ensemble 



Pa{L) = P{L) - Pbuik = - + lim J{L)/L (18) 

where J is the film grand potential per unit area. We can decompose the disjoining pressure into a contribution 
coming from the mean-field solution and a contribution coming from the field fluctuations which corresponds to the 
zero frequency van der Waals interaction 

pMF^^^_dJ^^ + lim J'^^{L)/L (19) 

Pf'^{L)=- ^ + lim r'''^{L)/L (20) 
After some straightforward but laborious algebra one finds that 



The mean-field value of the density within the film is given by 

Pmf(x) = 2pcosh(e/5V(x)) , (22) 
and in the linearized theory within the film this becomes 

Pmf(x) = 2p + pe2/32VXx)2 . (23) 
At the mid-plane z — L/2 of the film one has 

PM.(.= 2))-^-- = ^^^i^?(^- (24) 

The mid-plane pressure formula |^ for Poisson-Boltzmann theories with fixed surface charges or potentials, relates 
the disjoining pressure to the mid-plane mean-field density by 

P^'^ = kBT[pMF{z - ^) - pbulk] . (25) 

In fact, in theories of the type considered here with potentials at or near the film surface, one can show quite generally 
that the mid-plane formula holds generically as long as the external potential is zero in a finite interval containing 
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the mid-plane z = L/2. One can see directly, comparing Eq. ( |2l| ) and Eq. ( |25| ) with the linearized approximation 
Eq. (^), that the mid-plane formula is respected here. 

It is useful for what follows to work with the running variable a'{L) — D{L) tanh (■^)- The mean- field surface 
charge aMF{L) may be written in terms of a'{L) 

2pea'{L) 

crMF[L) = (26) 

m 

From Eq. ( p!?] ) we have 

a'iD (-oth (If))" (27) 

As i — !■ oo, a' takes its bulk surface value a* — limi^oo a'{L). As mentioned previously we are using a Gaussian 
field theory inside the film and so the mechanism generating the surface charge cannot be taken to be too strong, 
implying that a <C 1 which gives a* « a. As L decreases a'{L) decreases. However, since the expansion in Eq. ( [27| ) 
is in acoth(mL/2), the non-linear terms in this series must be taken into account when varying L. Indeed, for small 
L one can show that a'{L) ^ — mLln(i)/2, which tends to zero as L tends to zero. 
The mean- field disjoining pressure in terms of a'{L) is 



which, in terms of the mean-field surface charge reads 



2esinh2(li|^) ^ ' 

For the linearized Poisson-Boltzmann equation with constant surface charge ac one finds that p6| 

2 

Pmf = %:—r: ■ (30) 

2esinh^(i^i^) ^ ^ 

Hence, at the mean-field level, fitting the disjoining pressure at each value of L with an L dependent ac will reproduce 
the behavior of aMF{L) in the current theory. This result is also true for full non- linear mean- field theory and is not 
dependent on the quadratic field approximation in the bulk used here. 



IV. FLUCTUATION EFFECTS 



Evaluating the fluctuations about our mean-field solution yields 

^ ^rfxdy ^J'^ff^ l^^0'(x)0-(y)^-l / [3e{z){V4>'fd^-]- I (3em'cb'^ dK 



where 



-\pema\L) [ {S{z) + d{L ~ z))(j)'^dx , (31) 

a'{L) = D{L) tanh ( = ^^P +* e.p j-Pem) ^33) 
V 2 / 2p 



The main difference in Eq. ( p3| ) from the pure Gaussian theory of |23 is that the surface term a'{L) is now a function 
of the film thickness, whereas in the pure Gaussian theory it is a constant. The fluctuation term Eq. (|T]) may 
be evaluated using functional techniques ^ ^ ^ as it is a functional determinant, or by path integral 
techniques |^ . Using the results of |2^] one finds that the terms depending explicitly on a' and L that will contribute 
to the disjoining pressure and the surface charge are 



2 /. 

(3r'^^{L, a')^pj^XL + ^ j dkkXn [B{km)k + a\L) + 



, dfc fcln I 1 - ( B{km)k + a'{L)-VWTl_ ) (_2i^V^) 1 , (33) 
47r y 1 I B(fem)fc + a'(L)-f VFTT/ ' 
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where 



l + Aexp(-2p/i) 
= l-Acxp(-2p/.) ' 

with A = (eo — e)/((eo + e). The term /SJ^ui^jf in Eq. ( p3| ) is the van der Waals contribution to the the bulk grand 
potential per unit volume and is given by 

/34"uT -^Jdkk {V^l - k) . (35) 

The mean-field contribution to the bulk grand-potential per unit volume is simply PJ^^i^ = — 2/x since the mean- 
field solution in the bulk is just ip = 0; there are no surfaces to set up a mean-field potential. At this point we may not 
replace /i by p because we need to work to one-loop for consistency. The total bulk grand-potential per unit volume 
is thus given by 



/3Jbuik=-2Ai + ^ J dkk Vfc2 + 1 _ 



^ ^ 3 V 2to 

where A is a momentum space cut-off corresponding to a short distance cut-off a ~ 1/A and g = m^/Sirp ~ ttiIb- 
We see that indeed the expansion is in the dimcnsionlcss coupling g as asserted in section |l^ The bulk density of 
electrolyte is given by 

d TtV^ 1 

2/3 - /iT;-/3Jbuik = 2/^ -I- —Kvn? , 



which can be written as 



Zp, Z = 1-^ ( 1 + A 



From its definition in section |l| we then also have, to this order in g, that fi'^ = Z p'^. 

Substituting this result into Eq. (37) gives the well-known Debye expression for the bulk pressure 

/3Pbuik = Jbuik - 2p - — . (37) 

z47r 

In the field theoretic formulation used here the surface charge (on one surface) per unit area, cr, is given by 

a = e^;(exp(ie/3(^(z = 0)) , (38) 
which becomes, to the order of accuracy of the present treatment, 

a = e/4 exp {~ef3i^{z = 0)) (1 + iel3<j>'{z = 0) - ^e^ f3^(t)''^{z = 0)) , (39) 

where the average (•) in the above equation is over the fluctuations 0'. As the one-loop action in (/)' is quadratic, the 
average of the term linear in 0' in Eq. (|39|) is zero and we may then write 



2/iea'(L) ea'{L) df3J 



vdW 



|a'=o + 0(a'') , (40) 



m 2 da'{L) 



where we have used that a = mp*j^/2p = mp*^/2p and where we have kept only the leading order behavior of the 
surface charge in a' . The next order terms can be calculated and are finite though one needs to eliminate certain 
artificial divergences H. The formula Eq. ^ gives a surface charge susceptibility with respect to the conjugate 
variable a' . We note that the first term in on the right hand-side of Eq. (|4^) is simply the mean-field contribution to 
the surface charge (Jmf- 
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To 0{a') we obtain 



, , 2ea'(L) ( f kdk 



TO \ Sn J D^{k,m) 

+ ^ [kdk D-{k,m) \ ^ , 

16tt J D+{k, m) [D+{k, m)2 cxp{2LmVk^ + !) + £»_ (fc, to)2] / 

where D±{k,m) = kB{km) ± y/k^ + 1 . There is, however, a divergence in the first integral term in Eq. (^l|). This 
term corresponds to the van der Waals contribution to the charge of a bulk surface fibuik — limi^oo cr. This divergence 
can be regularized by choosing the fugacity fi to give the desired bulk electrolyte density p correct at one loop order. 
We define 

r = [ kdkf ^ ^ - 7^ , I (42) 

and then 

, (43) 

to \ OTT lOTT / 

where A is the ultra-violet cut-off introduced in the bulk calculation above. Note that for weak charging a* ^ a. 
Hence in terms of the physical variable p we obtain the divergence-free formula for Ubuik 

2ea* ( w? \ 

We remark here that in the case A 7^ 0, if the Stern layer thickness h is taken to zero then the term F in Eq. ( |4^ ) 
diverges. It is clear however that one cannot have a surface charge exactly at the interface between two media of 
different dielectric constants due to the presence of arbitrarily close image charges. The divergence in /i as /i — > in 
Eq. (Q) is thus a physical divergence and any model using a surface charge must place this surface charge away from 
a discontinuity in the dielectric constant. When A = no such divergence is present and B{km) ~ 1 leading to the 
simple formula 

CTbuik = (1 ^ I) ^mf{<x>) . (45) 
We notice that from Eq. ( ^7| ) that this equation may be written as 

/3-Pbuikea* 

o-buik = ■ (46) 

TO 

From Eq. (^J) we see that effect of electrostatic interactions is to reduce the surface charge from the value it would 
have had without interactions. This is because the excess anions left in the bulk pull the cationic surface charge into 
the bulk. The case where A > (i.e. eo > e) leads to B{km) > 1 for all k and, examining the integrand in the 
formula defining F, we find that F < and hence that positive A increases the surface charge above that of the case 
A = 0. This is to be expected physically as the image charges in this case attract the ions towards the medium of 
higher dielectric constant | p6[ . In the case where hm ^ 1 we may evaluate the integral defining F since the leading 
divergence as hm — + comes from the large k integration. We find 



g gA \ 



J f7A/i=^(oo) . (47) 



Again, we see that for A > (A < 0) the enhancement (reduction) of the surface charge which can be physically 
attributed to the presence of image charges. 

Finally the L dependence of the surface charge at 0(0;') is given by 



a{L) = aMF{L) f 1 - f - 5F + f / Wfc D^Km) 

\ 2 ^ 2j 7^+(fc,TO) [£)+(fc,TO)2exp(2LTOVFTT)-fZ?-(fc,TO)2] 



(48) 



where D±{k, to) = kB{km) ± vfc2+T 
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Thus, because a'{L) — > as L ^ 0, we find that ctmf, and hence the surface charge cr, vanishes as the film becomes 
thin. This is a physical result which is not picked up by a pure Gaussian theory . 

As previously stated the total disjoining pressure is composed of a mean-field contribution and a contribution 
coming from the fluctuations. We find the total of these two terms gives 



where 



^ 2esinh2(i^) I <7MF{L){l + a'{L)coth{^)) ^ ^ J ^ 'l^P{k)' ^ ' 



kB(km) + a'(L) - ^/k^ + 1 / ^ , , 

/ fc = ^ ^—^ cxp -Lmy/k^ + 1 . 50 

In the derivation of this formula the divergences which arise at intermediate stages of the calculation can, as in the 
case of the surface charge, be shown to cancel in the final result. 

The first term in Eq. ( |49| ) is simply the mean- field contribution to Pd, the second (in the same bracket) is the 
contribution coming from the dependence of J'"'^^ on L via the L dependence of a'. The last term is a form of 
screened zero frequency van der Waals contribution. In the absence of electrolyte we find that this last term gives 
a l/L-^ Casimir attraction across the film. In the presence of electrolyte this interaction is screened [Q and decays 
exponentially as exp(— 2mL) which is twice as quickly as the mean-field contribution to the disjoining pressure which 
decays as exp(— mL) for large L. However in the theory presented here, the fluctuation effects give an additional 
term which decays also as exp(— mL) and hence the fluctuations modify the long distance behavior of Pd- We note 
theories with fixed surface charge the long-range component of the disjoining pressure is not altered at two-loop level 
and this effect is specific to charge regularized models. The strength and sign of this long range modification of Pd 
is controlled the ratio ((t(L) — aMF{L)) / (Jmf{,L) which measures the relative deviation of the surface charge from 
its mean-field value. If the fluctuation effects enhance (reduce) the surface charge then the long range value of Pd is 
increased (decreased) from its mean-field value. The Casimir term is, however, is always attractive. 

It should be remarked at this stage that to this order, 0{g), there will be L-dependent contributions from the 
non-Gaussian interaction term in cj)'^ . These contributions can be shown to vanish as L ^ oo but their effect for finite 
L must be calculated. However, such a calculation requires the apparatus for the general perturbation theory to be 
developed which we shall present in a forthcoming paper [52| . 



V. DISCUSSION AND CONCLUSION 



We have systematically developed a theory for a thin film with a full non-linear surface charging mechanism while 
retaining the free field theory description for the bulk electrostatic fields. This corresponds to using linear Debye theory 
in the bulk with fugacity but with the fully interacting description of the sources for the charging mechanism of the 
surfaces. We have applied the theory to a model consisting of a triple layer system, shown in Fig. (|l|) , in which there 
is adsorption of cations on to the surface modelled by a surface fugacity ^\ and encoded in the dimensionless surface 
absorption strength parameter a = mp*^/2p where p and p*^ are the bulk density and cation density on the surface of 
a bulk region, respectively, and m is the Debye mass. At the surface there is a Stern layer of thickness h from which 
all ions are excluded. The film is of thickness L and the dielectric constants are e ^ 80eo in the film and eq outside 
the film. This model for a real surface is too simple but it encodes the important feature that the thermodynamic 
properties of the film are very strongly dependent on the detailed nature of the surface and its properties. This is due 
to two features: the charging mechanism which allows the surface charge to remain in equilibrium with the interior 
charges and the effect of image charges due to the discontinuity in the dielectric constant at the surface. First we 
determine the mean-field solution (pdz) using the non-linear surface operators as the source and then we use the 
Schrodinger kernel approach to calculate the partition function as an expansion in a and g — mis where Ib is the 
Bjerrum length. Much of the details of this approach have been discussed in an earlier paper [ p3[ . In this paper, we 
concentrate on the effects of the non-linear surface charging mechanism which leads us to introduce an effective, or 
running, surface charging parameter a'{L) and we analyze the behavior of the surface charge g(£) and the disjoining 
pressure Pd{L) on m, a and h. The formulas summarizing our findings are Eqs. ( ^ ) and (^). Examples of the 
solutions to these equations are shown in Figs. (|^) to (|l2|) and we now briefly discuss the salient features of these 
results. 

In Figs. (H) and (^) we show a'{L) as a function of L for various values of a and to, respectively. a'{L) controls 
the strength of the surface charging mechanism and from Eq. (|l^), and what follows, it is clear that a' [L) ^ as 
L — > which in turn causes umf and a to vanish also in this limit. 
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In Figs. (Q) to (^) we show the dependence of Pd on h,m,a, respectively. From all these figures we see that the 
characteristic collapse transition is evident but that its strength is very sensitive to the parameters. In particular, 
from Fig. we see that Pd decreases as h decreases, as we should expect since the image charges at the surface are 
repelling the cations and so reducing the surface charge. This effect can be seen directly in Fig.(||). The effect on 
Pd is due to the a dependent term in Eq. ( ^9| ) which arises because of the implicit L dependence of the free energy 
through its dependence on a'{L). Note that aMF{L) is independent of h as in this formulation the effect of image 
charges first comes in at the one- loop level. The dependence of Pd on m and a shown in Figs. (||) and (||) has the 
expected trends but again the height of the peak in Pd{L) is strongly dependent on the parameters which have been 
chosen to take values that can typically be achieved in experiment. In Fig. we show the mean-field component 
and one-loop contributions for typical parameter values. 

In Figs. (^) to ( pT| ) the behavior of the surface charge cr{L) is shown. We see from Fig (|9|) that a{L) decreases 
very strongly with L for small L and that as h decreases the effect strengthens and that for sufficiently small h the 
value of <j{L) has a zero and becomes negative. Ultimately, cr(L) — > as L — *■ and so must have a minimum value 
before turning towards the origin. These effects are due to the one-loop term in Eq. and can be seen clearly in 
Fig. ( |l2|) where the mean-field and one-loop contributions to ij{L) are separately shown. In Figs. (]l^) and ( |ll| ) the 
trends shown are as expected but, as in the case of Pd the magnitude of cr{L) is very sensitive to parameter values. 
The overall prediction is that (j{L) is strongly dependent on L and vanishes as L — > 0. 

In Fig. (H) we compare Pd for the linearized theory from (p3|) and the non- linear theory of this paper with 
h — 0.3 nm. Although the peak in Pd occurs in much the same place it is lower in the non-linear theory for this value 
of h. Since the peak height is strongly dependent on h we see that a quantitative prediction requires a realistic model 
for the surface. See Fig. (^. 

An important feature of these calculations is to note that the results are expressed as a series in both a or a'{L) and 
g ~ mis with partial resummations in some cases. The major approximation is to use the free field theory within the 
bulk. The object was to study the effects of the non-linear surface charging mechanism and we have shown that these 
effects are indeed strong and it is clear that any approach which omits them or assumes a constant surface charge will 
be incorrect. Some of the effects are strong and there are features, such as the minimum in a{L) and its change of 
sign, that must be studied further in the full non-linear theory to see if they are not artifacts of the approximation. A 
consistent control over spurious and artificial infinities must await the full perturbation theory. An example is eluded 
to in Eq. and what follows. We have indicated how to control such quantities here to the one-loop level and see 
that even here the analysis is rather delicate. There are, in principle, 0{g) terms from interactions within the bulk 
which will vanish as i — > oo but contribute finite L effects to Pd and cr but these terms are not expected to be large. 
We shall present an analysis of all these topics in a forthcoming paper ||3^ in which the full non-linear theory and its 
perturbation expansion will be studied. 
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FIG. 1: Cross section through a model soap film. The distance of closest approach of the ions to the surfaces (at = and 
^; = L) is the effective radius of the ionic species in solution. The dielectric constants and Debye masses as a function of the 
distance perpendicular to the film surface (ie.)as a function of the position on the z axis are also shown. 
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FIG. 2: q'(L) for Debye mass m = 0.2 nm ^ and h = 0.3 nm for values of a = 0.2,0.3,0.4,0.5 which controls the strength of 
the surface charging adsorption. 



12 



0.40 



0.30 



0.10 



0.00 





- - zsr -- - - '^ '^ '^ '^ 




/ ■ 




- 


/ / 






1 y 






/ / 
/ / 
/ ,•' / 

/ /'' / 
/ /■ / 




- 


■ ' ''' / 
' / 


m = 0.05 nm ' 






m = 0.1 nm"^ 






m = 0.2 nm"^ 




"7 
/ 







0.0 10.0 20.0 30.0 40.0 50.0 

L nm 



FIG. 3: ct{L) for a — 0.5 and Stem layer thickness h — 0.3 nm for values of Debye mass m — 0.05, 0.1, 0.2 nm ^ . 
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FIG. 4: The disjoining pressure given in Eq. (^) for Debye mass m = 0.2nm~^ and a = 0.5 for different values of the thickness 
of the Stern layer h = 0.1, 0.3, 0.5 nm. The sensitive dependence of Pd on h is evident as we should expect since the influence 
of the image charges increases rapidly as h decreases. 
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FIG. 5: The disjoining pressure given in Eq. 
mass m — 0.05, 0.1, 0.2 nm~^. 
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FIG. 6: The disjoining pressure given in Eq. ([lOl) for Debye mass m = 0.2 nm ^ and Stern layer thickness h = 0.3 nm for 
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FIG. 7: The disjoining pressure given in Eq. ( p^ showing the mean-field contribution and the one-loop (0(<?)) contributions 
as a function of Stern layer thickness h = 0.3. 
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FIG. 8: The disjoining pressure given in Eq. for m = 0.2 nm~^ and a = 0.5 for the linearized theory from (|2^) and 

the non-linear theory of this paper for which h = 0.3 nm. Although the peak in Pd occurs for much the same value of L it is 
lower in the non-linear theory for this value of h. Since the peak height is strongly dependent on h we see that a quantitative 
prediction requires a realistic model for the surface. See Fig. (^. 
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FIG. 9: The surface charge a in millicoulombs in Eq. (^) for Debye mass m — 0.2nm~^ and a — 0.5 for different values of 
the thickness of the Stern layer h — 0.1, 0.3, 0.5 nm. We see that a decreases with L and this effect is enhanced as h becomes 
smaller as we should expect since the image charges have greater influence. However, even though a is constrained to vanish 
we see that for small e noug h h it changes sign and so must have minimum at very small L before turning towards zero. This 
effect is evident in Fig (|l2[). Whether or not this effect is an artifact will be a subject of further study. 
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FIG. 10; The surface charge a in millicoulombs in Eq. 
Debye mass m = 0.05, 0.1, 0.2 nm^^. 



for a = 0.5 and Stern layer thickness h = 0.3 nm for values of the 
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FIG. 11: The surface charge a in miUicoulombs in Eq. Mm for Debye mass m = 0.2 nm ^ and Stern layer thickness h — 0.3 nm 
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FIG. 12: The surface charge a in miUicoulombs in Eq. (|4^) showing the mean-field contribution and the one-loop (0(gi)) 
contributions as a function of Stern layer thickness h — 0.1,0.3,0.5 nm. The one-loop contribution is negative and has a 
minimum before turning to zero as it must since a'{L) vanishes as L — > 0. 
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